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Abstract. A study of the angular distributions of leptons from decays of J ftp's produced in p-C and p-W 
collisions at ^/s = 41.6 GeV has been performed in the J/ip Feynman-a; region —0.34 < xf < 0.14 and 
for J/ip transverse momenta up to 5.4 GeV/c. The data were collected by the HERA-B experiment at 
the HERA proton ring of the DESY laboratory. The results, based on a clean selection of 2.3 • 10® J/ip's 
reconstructed in both the e^e" and ijT decay channels, indicate that J/ip's are produced polarized. 
The magnitude of the effect is maximal at low pr- For pr > 1 GeV/c a significant dependence on the 
reference frame is found: the polar anisotropy is more pronounced in the Collins-Soper frame and almost 
vanishes in the helicity frame, where, instead, a significant azimuthal anisotropy arises. 
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1 Introduction 

This paper presents a new measurement of the angu- 
lar distribution of leptons from J/^'s produced inclu- 
sively in proton-nucleus collisions at centre-of-mass en- 
ergy s/s — 41.6 GcV. The data were collected by the 
DESY experiment HERA-B and covered the kinematic 
ranges —0.34 < xp < 0.14 in the Feynman-x variable and 
< < 5.4 GeV/c in transverse momentum. In this do- 
main, the average fraction of J mesons coming from Xc 
and %p' decays has been determined as ~ 27% Most 
previous analyses were based on the choice of one specific 
definition of the polarization frame and were limited to the 



measurement of the polar angle distribution, from which 
the so-called "polarization" parameter is extracted. The 
present measurement includes for the first time a system- 
atic comparison of the results obtained for the full decay 
angular distribution in three different reference frames - 
and significant differences are found between them. Some 
of the results are presented separately for the two target 
materials used in the experiment (carbon, A — 12, and 
tungsten, A = 184), leaving open the possibility that the 
nuclear medium may affect the observed decay kinemat- 
ics (for example as a consequence of a varying mixture 
of J/ip^s from decays of heavier charmonium states and 
direct J/ip's). The analysis is based on almost the same 
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sample used in the measurement of the J /tp kinematic dis- 
tributions described in [3] , with a total of about 83000 and 
143000 J/ip's reconstructed, respectively, in the dimuon 
and dielectron decay channels (excluding only the small 
fraction of data collected with titanium). The reader is 
referred to that paper for a description of detector, data 
taking, trigger, selection criteria and Monte Carlo simula- 
tion. 

The next section explains the definitions and con- 
ventions used in the measurement (Sect. [2|). The results 
are presented in Sect. [3] and discussed in the conclusions 
(Sect. SI). 



2 Definitions 

In any chosen reference system, the most general (parity 
conserving) form of the two-lepton decay angular distri- 
bution of a J/f/; is 



and 



dN 



d{coa6) d(j) 



(X 1 + Xgcos^e + X4,sm^ 9coa{2(f>) 
+ \ed> sin(20) coach, 



(1) 

where 9 is the angle between the direction of the posi- 
tive lepton and a chosen polarization axis in the J /tp rest 
frame. </>, the corresponding azimuthal angle, is defined 
with respect to the plane of the colliding hadrons; Ag, Ag^ 
and A^ are the quantities to be measured^ The param- 
eter Xg is usually called "polarization" ; more generally, 
any observed deviation of at least one of the three param- 
eters from zero would be the indication of polarized J/t/i 
production. Non-zero values of A^ and Xq^ indicate an az- 
imuthal anisotropy of the distribution, which may arise as 
a consequence of specific choices of the reference frame. 

The following different definitions of the polarization 
axis are commonly used in the literature: 

— the direction of the beam momentum (as seen in the 
J/ip rest frame) in the so-called Gottfried- Jackson 
(GJ) frame @]; 

— the direction of the J /ip in the center-of-mass system 
of the colliding hadrons (helicity frame, HX); 

— the bisector between the directions of the beam mo- 
mentum and of the opposite of the target momentum 
(as seen in the J/ip rest frame) in the CoUins-Soper 
(CS) frame [5]. 

The results of the present analysis are presented in all of 
these frames. 

The values of Xg, X^ and Xg^ are extracted from the 
data by considering single-variable projected angular dis- 
tributions. When averaged over (p and cos 6*, respectively, 
the angular distribution of the decay leptons is reduced to 
the forms 

dN 



c?(cos 9) 



cx 1 + Ae cos 9 



(2) 



^ The sign of the measured Xg^ depends on the chosen orien- 
tation of the y axis (the one perpendicular to the plane of the 
colliding hadrons). In the convention adopted for the present 
analysis, the y axis is oriented as the vector product of the 
beam momentum and the target momentum, pbcam x ptargct. 



dN , 2A<A 



(3) 



One possible way of determining the value of Xg^p is to 
define the variable 6g as 



jTT for cos 6* < 
- 2: for cos 61 > 



and measure the distribution 



dN ^ , V2Xg^ 

-TT- « 1+ q , , COS' 

dchg o -f Ag 



(4) 



(5) 



3 Results 

The efficiency-corrected single-variable angular distribu- 
tions averaged over the accepted phase space (—0.34 < 
xf < 0.14, < pt < 5.4 GeV/c) are shown in Fig. [1] with 
statistical uncertainties only. The results are given sepa- 
rately in the two decay channels (/i"*"/!^, e^e~) for car- 
bon and tungsten target data combined. As can be seen, 
the measured distributions follow the correct symmetric 
and/or periodic behaviour expressed by the formulas in 
Eqs.[21[Slaiid[ni indicating a good level of reliability of the 
acceptance and efficiency correction procedures; the chi- 
square probabilities obtained by fitting the distributions 
are on average 15% and 50% in the muon and electron 
channels, respectively. The values of the output parame- 
ters of the fits are listed in Table[TJ The difference between 
the values measured in the two channels is always less than 
1.3(7. All results presented hereafter are averages of muon 
and electron measurements. We remark that a prelimi- 
nary value of —0.35 ± 0.04 was used for the effective Xg 
(CS frame) in the evaluation of systematic uncertainties 
of a previous HERA-B analysis [2J; such value deviates 
slightly, but well within errors, from the combined value 
in Table [H 

The final results including the estimated systematic 
uncertainties are displayed in Fig. [2] as a function of the 
transverse momentum and Feynman-x of the J/V'- As be- 
fore, the two target data samples have been combined. 
The corresponding numerical values are listed in Tables [21 
OandUl where (pr) and (xp) indicate averages over the 
J/ip's reconstructed in a given bin. The bin boundaries are 
defined by the following lists: 0, 0.5, 1.0, 1.5, 2.5, 5.4 GeV/c 
(0,0.2,0.4,0.6,0.8,1.0,1.5,2.5,5.4 GeV/c for the helicity 
frame) for pr and -0.34, -0.15, -0.10, -0.05, 0, 0.05, 0.14 
for X F ■ The systematic errors have been evaluated with the 
procedure already described in our report on the measure- 
ment of the J/ip kinematic distributions [3], taking into 
account the impact of signal selection and optimization, 
signal counting method, differences in acquisition condi- 
tions, and the kinematics of the MC generation. Addi- 
tional systematic tests consisted in fitting the angular dis- 
tributions excluding the angular ranges with the lowest 
efficiency. Statistical and systematic errors are obviously 
correlated from frame to frame. The systematic errors in 
different px/xp bins are partly correlated. 
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Fig. 1. EfRciency-corrected distributions of the angular variables cos6' (a, b), (ji (c, d) and (j}e (e, f) measured in the muon (left 
column) and electron (right column) decay channels of the J/tp over the whole visible phase space. The results obtained in the 
CoUins-Soper, Gottfried-Jackson and helicity frames are represented, respectively, by black circles, white squares and asterisks. 
The errors are only statistical. For a better visualizations, the CoUins-Soper and helicity distributions are displayed with their 
actual values shifted by a constant. The distributions are fitted with the curves of Eqs. |2l [3] and [5] 



Table 1. Output parameters obtained by fitting the distributions shown in Fig.[T]with the curves of Eqs. [2] [3] and [S] The errors 
in the parameters reflect only statistical uncertainties of the distributions. The errors are correlated from one frame to another. 
The systematic uncertainties for the average values are in all frames of the order of 0.05 for Ae, 0.02 for and 0.015 for Ae0. 



Frame /channel 


As 


A^ 


Ae0 


CS/e+e- 
CS/avg. 


-0.296 ±0.029 
-0.383 ±0.061 
-0.313 ±0.026 


-0.0194 ±0.0051 

-0.022 ±0.011 
-0.0199 ±0.0046 


0.0158 ±0.0049 
0.0195 ±0.0096 
0.0168 ± 0.0043 


GJ/e+e" 
GJ/avg. 


-0.185 ±0.021 
-0.256 ±0.051 
-0.195 ±0.019 


-0.0400 ± 0.0051 

-0.031 ±0.011 
-0.0385 ± 0.0046 


0.0433 ± 0.0051 
0.058 ±0.010 
0.0463 ± 0.0045 


RX/fj.+fi- 
nX/e+e- 
HX/avg. 


-0.115 ±0.012 
-0.092 ±0.027 
-0.111 ±0.011 


-0.0714 ±0.0055 

-0.075 ±0.012 
-0.0720 ± 0.0050 


0.0049 ± 0.0049 
0.0161 ±0.0094 
0.0073 ± 0.0043 
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Fig. 2. The parameters Xe (a, b), (c, d) and Ae^ (e, f) measured as functions of the average reconstructed pr (left column) and 
xf (right column). The results obtained in the CoUins-Sopcr, Gottfried- Jackson and helicity frames are represented, respectively, 
by black circles, white squares and asterisks. The vertical errors bars represent quadratic sums of statistical and systematic 
uncertainties. The horizontal bars indicate the adopted binning. 
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The results indicate an anisotropy of the J /if) decay an- 
gular distribution, visible in either its polar or azimuthal 
projections (in the CS and HX frames, respectively), or 
in both (GJ frame). Moreover, there is a definite hierar- 
chy in the magnitudes of the parameters Xq and A^. In 
particular, the polar anisotropy (Ag < 0) increases when 
going from the HX to the CS frame, while the azimuthal 
parameter A^ changes following a reversed order. Both 
parameters have in-between magnitudes in the GJ frame. 
A kinematic dependence characterizes the results. For ex- 
ample, the magnitude of the polarization parameter \g in- 
creases with decreasing pT ■ This low-px effect is the same 
in the three frames - as expected from the fact that all 
frames coincide at px — - whereas different polariza- 
tion magnitudes are measured at higher px^ following the 
already mentioned hierarchy. 

By performing a target-dependent analysis, gener- 
ally small differences between the polarization parame- 
ters measured with carbon and tungsten data have been 
found and taken into account in the systematic errors of 
the already shown combined measurements. However, the 
results for the parameter \g show a slightly significant 
target dependence in the lowest pT bin. The polarization 
measured in the CS frame as a function of px and Xi? is 
shown in Fig. [3] for the two target materials; the num- 
bers are listed in Table O In the first px bin, carbon and 
tungsten results differ by about 3cr at the statistical level; 
the significance of the difference is reduced to the 2a level 
when also the systematic errors are included in the com- 
parison. 

Since the largest polarization effects are seen at low 
px^ special attention has been devoted to the investiga- 
tion of the systematic variations in the value of Xq for 
< pt < 0.5 GeV/c. In fact, the determination of the 
low-pT efficiency is especially sensitive to the Monte Carlo 
description of the detector region close to the beam pipe. 
Various tests have been performed, which consisted in re- 
ducing the geometrical acceptance around the beam pipe, 
changing the beam direction in the Monte Carlo and test- 
ing different selection criteria for the low momentum lep- 
tons. As an extreme check, the acceptance determination 
has been varied to reproduce the detector conditions of 
acquisition periods not correlated to the considered sam- 
ple. Final systematic uncertainties of ~ 0.1 and ~ 0.2 have 
been evaluated for Xq for the carbon and tungsten targets, 
respectively (Table [5]). 



4 Conclusions 

HERA-B has measured the two-lepton decay angular dis- 
tribution of J/ip^s produced inclusively in proton- nucleus 
collisions, using the decay channels e~^e~ and ^~^fj,~ . The 
distributions of the polar and azimuthal angles have been 
determined in three different polarization frames. The re- 
sults can be summarized as follows. 

— In the observed phase space, Xg is negative, indicat- 
ing that the J/ip^s are produced with a preferred spin 
component along the reference axis. 



— There is a definite hierarchy for the values of the decay 
angular parameters measured in different frames: the 
polar and azimuthal parameters satisfy the relations 

|Ae(HX)| < |Ae(GJ)| < |Ae(CS)| 
|A4HX)| > |A4GJ)| > |A4CS)|, 

while Xg^ is significantly different from zero only in the 
GJ frame. 

— The polarization effects depend on the kinematics of 
the J/ip. In particular, the polar anisotropy increases 
with decreasing and is maximal in the limit px 0. 
The different results obtained in the three frames - 

in terms of both polar and azimuthal distributions - are 
an example which shows that an analysis limited to only 
one frame and one polarization parameter is in general in- 
complete. For example, the present measurement of only 
the polarization parameter Xg (i.e. ignoring A^) in the HX 
frame for px > ^ GeV/c may be misunderstood as a sig- 
nificant indication of unpolarized J/ip production. 

Among existing measurements of the parameter Xg, 
E866 [6J (p-Cu at = 38.8 GeV) has measured in the CS 
frame a pT-iudependent polarization consistent with zero, 
while the CDF Run II [7] data {pp at = 1.96 TeV) indi- 
cate a negative polarization in the HX frame increasing in 
magnitude with increasing p^- These results have been ob- 
tained in kinematic ranges (E866: xp > 0.25, CDF: px > 5 
GeV/c) which have no overlap with the HERA-B data and 
between each other. The three results are therefore not in 
contradiction; their comparison has to be interpreted as 
a further indication that the observed polarization effects 
change with varying kinematic conditions of the produced 
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Table 2. Values of the parameters \e, X<p and Xg^ measured in the CoUins-Soper frame as functions of the average reconstructed 
Pt and xf for combined carbon and tungsten data. The errors are statistical and systematic. 



(Pt) (GeV/c) 




X^ 


Xe4, 


0.35 
0.79 
1.25 
1.89 


-0.56 ±0.07 ±0.21 
-0.28 ± 0.05 ± 0.13 

-0.199 ±0.046 ±0.083 
-0.164 ± 0.049 ±0.039 
-O.l.")!) = O.OlO = O.O.'il 


-0.026 ±0.013 ±0.027 
-0.017 ±0.010 ±0.020 

-0.031 ±0.012 ±0.018 
-0.028 ±0.017 ±0.023 

-().()():! = ().0:!o = D.OGl 


0.042 ± 0.020 ± 0.022 
-0.001 ± 0.014 ± 0.022 

0.016 ±0.014 ±0.021 
0.046 ± 0.015 ± 0.017 
0.012 = 0.1)21 = I). Oil 




{xf) 


Xe 




Ae0 


-0.193 

-0.123 
-0.074 
-0.026 
0.021 
0.067 


-0.202 ± 0.072 ± 0.116 

-0.179 ±0.050 ±0.068 
-0.296 ± 0.052 ± 0.068 
-0.284 ± 0.051 ± 0.094 
-0.10 ±0.11 ±0.15 
-0.51 ±0.27 ±0.22 


-0.014 ± 0.019 ± 0.032 
-0.017 ±0.014 ±0.014 
-0.024 ±0.010 ±0.012 
-0.025 ±0.011 ±0.017 
-0.072 ±0.013 ±0.029 
-0.030 ± 0.025 ± 0.047 


0.039 ± 0.023 ± 0.057 

-0.003 ±0.019 ±0.024 
0.009 ±0.016 ±0.018 
0.021 ±0.013 ±0.025 
0.040 ± 0.018 ± 0.044 
0.007 ± 0.039 ± 0.076 



Table 3. Values of the parameters Xg, X^ and Xg^ measured in the Gottfried- Jackson frame as functions of the average 
reconstructed pr and xf for combined carbon and tungsten data. The errors are statistical and systematic. 



(pt) (GeV/c) 


Xg 


A^ 


Xg^ 


0.35 
0.79 
1.25 
1.89 
3.11 


-0.53 ±0.06 ±0.19 
-0.219 ±0.033 ±0.091 
-0.096 ±0.033 ±0.059 
-0.075 ± 0.025 ± 0.031 
-0.018 ±0.041 ±0.033 


-0.029 ±0.014 ±0.028 
-0.022 ±0.010 ±0.021 
-0.048 ±0.014 ±0.018 
-0.064 ± 0.017 ± 0.022 
-0.171 ± 0.036 ± 0.054 


0.069 ± 0.020 ± 0.020 
0.046 ±0.014 ±0.027 
0.049 ± 0.017 ± 0.032 
0.058 ± 0.018 ± 0.037 
0.063 ± 0.026 ± 0.038 




{xf} 


Xe 


A^ 


Xeci> 


-0.193 
-0.123 
-0.074 
-0.026 
0.021 
0.067 


-0.096 ± 0.039 ± 0.075 
-0.120 ±0.033 ±0.063 
-0.151 ± 0.028 ± 0.063 
-0.141 ±0.031 ±0.070 
-0.140 ±0.057 ±0.084 
-0.10 ±0.10 ±0.10 


-0.023 ±0.021 ±0.026 
-0.065 ±0.015 ±0.020 
-0.050 ± 0.012 ± 0.020 
-0.041 ±0.011 ±0.023 
-0.084 ±0.013 ±0.029 
-0.014 ± 0.026 ± 0.039 


0.050 ± 0.024 ± 0.063 
0.046 ± 0.020 ± 0.031 
0.059 ± 0.018 ± 0.024 
0.083 ±0.014 ±0.029 
0.088 ±0.019 ±0.045 
0.014 ± 0.049 ± 0.072 



Table 4. Values of the parameters Xg, X^, and Ae^ measured in the helicity frame as functions of the average reconstructed pt 
and Xf for combined carbon and tungsten data. The errors are statistical and systematic. 



{pt) (GeV/c) 


Xe 


A,^ 


Xe<j> 


0.14 
0.31 
0.51 
0.71 
0.90 
1.25 
1.89 
3.11 


-0.49 ±0.07 ±0.20 
-0.31 ±0.05 ±0.11 
-0.153 ±0.031 ±0.096 

-0.164 ± 0.026 ±0.085 
-0.093 ± 0.028 ±0.075 
-0.017 ± 0.020 ±0.059 
-0.037 ± 0.021 ± 0.037 
-0.001 ±0.037 ±0.016 


-0.045 ± 0.029 ± 0.047 
-0.005 ± 0.021 ± 0.041 
-0.025 ± 0.017 ± 0.034 

-0.058 ±0.018 ±0.029 
-0.031 ±0.017 ±0.025 
-0.091 ±0.022 ±0.021 
-0.090 ± 0.026 ± 0.023 
-0.243 ± 0.055 ± 0.063 


0.075 ± 0.043 ± 0.026 
0.075 ± 0.028 ± 0.027 
0.022 ± 0.026 ± 0.028 

-0.021 ± 0.023 ± 0.030 
0.051 ±0.023 ±0.031 

-0.012 ±0.017 ±0.034 
0.014 ± 0.017 ± 0.042 
0.032 ± 0.027 ± 0.062 




{xf) 


Xg 


A<^ 


Xg^f, 


-0.193 
-0.123 
-0.074 
-0.026 
0.021 
0.067 


-0.092 ± 0.037 ± 0.067 
-0.117 ± 0.026 ± 0.050 
-0.070 ± 0.024 ± 0.043 
0.014 ± 0.019 ± 0.041 

-0.075 ± 0.023 ± 0.043 
-0.064 ± 0.041 ± 0.049 


-0.037 ±0.023 ±0.025 
-0.067 ± 0.017 ± 0.027 
-0.085 ± 0.015 ± 0.030 
-0.112 ±0.023 ±0.034 

-0.100 ±0.027 ±0.040 
-0.028 ± 0.029 ± 0.047 


0.059 ± 0.025 ± 0.062 
0.024 ± 0.021 ± 0.035 
0.058 ± 0.016 ± 0.029 
0.073 ± 0.016 ± 0.033 

0.042 ± 0.020 ± 0.048 
-0.085 ± 0.044 ± 0.073 
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Table 5. Values of the parameter \g measured in the CoUins-Soper frame as a function of the average reconstructed pr and 
xf for carbon (C) and tungsten (W) data. The errors are statistical and systematic. 



(Pt> (GeV/c) \g (C) \e (W) 



0.35 -0.38 ±0.10 ±0.11 -0.81 ±0.10 ±0.20 

0.79 -0.271 ± 0.062 ± 0.078 -0.32 ± 0.09 ± 0.17 

1.25 -0.192 ±0.054 ±0.066 -0.22 ± 0.09 ± 0.13 

1.89 -0.147 ±0.062 ± 0.052 -0.227 ± 0.081 ± 0.076 

-O.Kii = I). 01!) = ().():!<) -O.iUi = 0.1)72 = I). 038 



{xf) 


\e (C) 


\e (W) 


-0.193 


-0.097 ± 0.090 ± 0.099 


-0.32 ± 0.12 ± 0.08 


-0.123 


-0.176 ±0.054 ±0.062 


-0.18 ±0.15 ±0.10 


-0.074 


-0.304 ± 0.079 ± 0.067 


-0.29 ±0.07 ±0.10 


-0.026 


-0.236 ± 0.062 ± 0.097 


-0.37 ±0.09 ±0.10 


0.021 


0.03 ±0.15 ±0.15 


-0.21 ±0.17 ±0.10 


0.067 


-0.71 ± 0.31 ± 0.23 


-0.02 ± 0.52 ± 0.10 



0.2 - 

-- 

-0.2 - 
'< -0.4 — 
-0.6 - 

-0.8 — 
-1, 



0.5 1 



1.5 



2 2.5 3 
Pj [GeV/c] 



3.5 



(a) 

4 ' 4!5 



0.2 - 


-0.2 - 



-0.4 - 
-0.6 - 
-0.8 - 

"io.3 -0.25 -0.2 -0.15 -0.1 -0.05 



(b) 

0.05 0.1 



Fig. 3. The parameter Xe measured in the CoUins-Soper frame in carbon (black points) and tungsten (white points) as functions 
of the average reconstructed pr (a) and xf (b). The vertical error bars represent quadratic sums of statistical and systematic 
uncertainties. The horizontal bars indicate the adopted binning. 



